Preparation of atomic Fock states by trap reduction 
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We describe the preparation of atom-number states with strongly interacting bosons in one di- 
mension, or spin-polarized fermions. The procedure is based on a combination of weakening and 
squeezing of the trapping potential. For the resulting state, the full atom number distribution is 
obtained. Starting with an unknown number of particles Ni, we optimize the sudden change in 
the trapping potential which leads to the Fock state of Nf particles in the final trap. Non-zero 
temperature effects as well as different smooth trapping potentials are analyzed. A simple criterion 
is provided to ensure the robust preparation of the Fock state for physically realistic traps. 

PACS numbers: 32.80.Pj, 05.30.Jp, 05.30.Fk, 03.75.Kk 



I. INTRODUCTION 



Given the importance of photon statistics in quantum 
optics, the field of atom statistics is expected to develop 
vigorously in atom optics, fueled by the current ability 
to measure the number of trapped ultracold atoms with 
nearly single-atom resolution and without ensemble aver- 
aging of fluctuations [E[3,S| . Among the possible atomic 
distributions, pure atom number (Fock) states form a 
fundamental basis and hold unique and simple properties 
that make them ideal for studies of quantum dynamics 
of few-body interacting systems [3], precision measure- 
ments or quantum information processing [H, 0, H|- 
Efficient and robust creation, detection and manipula- 
tion of atom number states are thus important goals in 
atomic physics. Several approaches have been proposed 
and explored recently with theoretical and experimental 
work leading to sub-Poissonian and, in the limit, number 
states, such as atomic tweezers [l(|, interferometric 
methods fl . |Tll| . Mott insulator states [13], or atomic 
culling [U ll3j]. None of these methods is so far fully sat- 
isfying if the individual atoms have to be addressed (a 
problem of the Mott insulator states in optical lattices), 
and if an arbitrary number of atoms is to be produced 
reliably and with small enough variance for the trapped 
atom number, so further research is still required. 

In a previous paper [l4| a method was proposed in 
which the trapping potential is simultaneously weakened 
and squeezed so that the final trap holds a desired num- 
ber state. For a Tonks-Girardeau gas, it was shown that 
this mixed trap reduction yields optimal results, even 
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when the process is sudden. From the expression of the 
number variance, a simple criterion for optimal perfor- 
mance was obtained, namely, that the subspace spanned 
by the occupied levels in the initial trap configuration 
contains the subspace of the bound levels in the final trap. 
Starting from an unknown number of particles trapped 
at zero temperature, the mixed trap reduction method 
assures that the final state indeed corresponds to the 
desired Fock state by avoiding the momentum or posi- 
tion space truncations inherent in pure squeezing or pure 
weakening (the latter being called "culling" in [131]). 
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141 ] the potential traps considered for simplic- 
ity were finite square wells, so doubt could be cast on the 
validity of the results in actual smooth traps. Other limi- 
tations were the consideration of zero temperature initial 
states, and a statistical analysis limited to the first and 
second moment of the number distribution. In this paper 
we overcome these shortcomings by studying the mixed 
trap reduction process using smooth potentials, states 
with finite temperature, and the full number distribu- 
tion. 



II. THE TONKS REGIME AND POLARIZED 
FERMIONS 

The strongly interacting regime of ultracold bosonic 
atoms can be described by the so-called Tonks-Girardeau 
(TG) gas [HI, which is achieved at low densities and/or 
large one-dimensional scattering length 0, It has 
been argued [l4[ that this regime is optimal for the cre- 
ation of atomic Fock states by mixed trap reduction. 

The TG gas and its "dual" system of spin-polarized 
ideal fermions behave similarly, and share the same 
one-particle spatial density as well as any other local- 
correlation function, while differ on the non-local corre- 
lations. 

The fermionic many-body ground state wavefunc- 
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tion of the dual system is built at time t = as a 
Slater determinant for N particles, ipp (x\ , . . . , x^ i ) = 
7^T det ^fc=i^™(a;fc) ; where ^(a;) is the ra-th eigen- 
state of the initial trap, whose time evolution will 
be denoted by ip n (x,t) when the external trap 
is modified. The bosonic wave function, sym- 
metric under permutation of particles, is obtained 
from ipp by the Fermi-Bose (FB) mapping [TH, [l8[ 
tp(xi,...,x Ni ) = A(xi, . . . ,x Ni )ip F (xi, . . . ,x Ni ), where 
A = J[ 1< j <k<Ni sgn(x fc - Xj) is the "antisymmet- 
ric unit function". Noting that \tp\ 2 = \ip F \ 2 it is 
clear that both systems obey the same counting statis- 
tics. Moreover, since A does not include time explic- 
itly, the mapping is also valid when the trap Hamilto- 
nian is modified, and the time-dependent density pro- 
file resulting from this change can be calculated as 
[H[ p(x,t) = Njf \tp(x,x 2 , ■ ■ ■ ,x Nz ;t)\ 2 dx 2 • ■ -dx Ni = 
J2n=i \ fn(x,t)\ 2 . By reducing the trap capacity (maxi- 
mum number of bound states and thus particles that it 
can hold in the TG regime) some of the Ni atoms initially 
confined may escape and only TV will remain trapped. 

To determine whether or not sub-Poissonian statistics 
or a Fock state are achieved in the reduced trap we need 
to calculate the atom- number fluctuations. 



III. THE SUDDEN APPROXIMATION: FULL 
COUNTING STATISTICS 

We shall now describe the preparation of Fock states 
by an abrupt change of the trap potential to reduce its 
capacity. Consider a trap with an unknown number of 
particles Ni, which supports a maximum of Cj bound 
states. Generally Ni is smaller than the capacity of the 
trap Ci . The trapping potential is abruptly modified to a 
final configuration of smaller capacity Cf. Similarly the 
final number of trapped particles will be Nf < Cf. We 
are interested in the optimal potential change such that 
Nf = Cf to prepare the atomic Fock state \Nf — Cf). 

Let a — i, f stand for initial and final configuration. 
The Hilbert space associated with the Hamiltonian of a 
particle moving in any realistic trap V a , is the direct sum 
H a = B a © Ha of the subspace spanned by the bound 
states B a — {\<p")\j = 1, ■ ■ ■ , C a }, and that of scattering 
states lZ a = {\Xk)\k e Consider the projector onto 
the final bound states, Bf defined as 

Cf 

Within the TG regime and for spin-polarized fermions, 
the asymptotic mean number and variance of trapped 
atoms are [l4{ 

(N f ) = Tr(AiAf) (2) 

and 

a% f = TrfAiAf-^Af) 2 ], (3) 



where 

Ni 

A i = ^K)R| (4) 

71=1 

is the projector onto the bound subspace occupied by 
the initial state. We may thus conclude that trap reduc- 
tion can actually lead to the creation of Fock states with 
{Nf) = Cf and af^ f — quite simply when the initial 
states span the final ones, 

A f c A,. (5) 

In fact the full atom number distribution [2(| is accessible 
in the atom culling experiments [l[ and we next focus our 
attention on it. Consider the characteristic function of 
the number of particles in the bound subspace of the final 
trap, 

F(9) = Tr[pe wl(&Kf ]. (6) 

Following [2l[, the atom number distribution can be ob- 
tained as its Fourier transform, 

P(n) = ^f e- inB F{6)d6, (7) 

with n = 1, . . . , Ct. 

The characteristic function of spin-polarized fermions 
or a Tonks-Girardeau gas restricted to a given subspace 
was studied in [13, HJH3l ■ Using the projector for the 
bound subspace in the final configuration, 

F(6>)=det[l + (e 4e -l)A / A, ; ]. (8) 

For computational purposes it is convenient to use the 
basis of single-particle eigenstates |<p£i) which spans the 
final bound subspace so that F(ff) = detA, where A is a 
Cf x Cf matrix with elements 

A nm = S nm + [exp(i6) - l](^|Ai|^4>- (9) 

Clearly, if A/ c A i; A nm = exp(i6)6 nm , F(9) = 
exp(i8Cf) and the Kronecker-delta atom-number distri- 
bution associated with the Fock state \N — Cf) is ob- 
tained, 

p(n) = 5 nt c r (10) 

For completeness we note that the cumulant-generating 
function log F(9) admits the expansion 

logF(0) = f> n ^f , (11) 

n—l 

from which the mean k\ = (Nf) in Eq. @ and variance 
K 2 = 0jv, in Eq. (J3j) are just the first two orders. 

Different regimes of interactions for ultracold Bose 
gases in tight waveguides can be characterized by a sin- 
gle parameter 7 = mgiuL/h 2 N, where giu is the one- 
dimensional coupling strength, L the size of the system, 
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FIG. 1: Schematic potential change for trap weakening (a) 
and squeezing (b), relative to the initial width Li and depth 
Vi, for a = 0.05L. 



m and N the mass and number of atoms respectively. 7 
can be varied [17| allowing to explore the physics from 
the mean- field regime (7 <C 1) to the TG regime (7 3> 1) 
[l6| . We note that for the system to remain in the TG 
regime, it suffices to keep or decrease the density, since 



7/ = li 



(12) 



In particular, trap weakening clearly leads to a reduction 
of the density so that the system goes deeper into the 
TG regime. 



IV. DEPENDENCE ON THE TRAPPING 
POTENTIAL 

In this section we shall discuss the efficiency of the 
trap-reduction procedure at zero-temperature, focusing 
on the relevance of the shape of the confining potential. 
In particular, instead of the idealized square potentials 
used in [l3|, [3| we shall study here the family of "bath- 
tub" potentials 



V a {x', V a , L a , a a ) — ~2^ a 



1 — tanh 



\x\ - L a 



(13) 



as well as the inverted Gaussian potential 

V(ar) = -V a e~^. (14) 

For the bathtub L a and V a play respectively the role of 
the width and depth of the trap, while a a is an addi- 
tional parameter describing the smoothness of the po- 
tential trap. 

The spectrum and eigenfunctions can be found nu- 
merically by a standard technique, first differencing the 
Hamiltonian and then diagonalizing the tridiagonal ma- 
trix obtained by such difference scheme [25j . 

For the bathtub potential, a given U a — 2mV a L a /h 2 
and a a = <J a /L a defines a family of isospectral poten- 
tials. In dimensionless units, their eigenvalues are the 
same and so are their eigenfunctions. In the limit of a 
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FIG. 2: Atom number distribution p(n) in a trap reduction 
setup combining both squeezing and weakening techniques, 
for a bathtub potential given by Eq, (ll3|l with a a — 0.03L a 
and the parameters of the trap being iV, = d — 100 and 
C f = 10 (Ui = 10 4 tt 2 , Uf = 10 2 tt 2 ). Plot (a) corresponds 
to almost pure squeezing, with Lf/Li = 0.04, while plot (c) 
corresponds to pure weakening, Lf/Li — 1, showing that it 
is the combination of both techniques, Lf/Li = 0.5, plot (b), 
the most efficient way for our purposes. 



square potential U a « C a 2 n 2 . For the Gaussian poten- 
tial a single parameter U a = VaS^ defines an isospectral 
family. 

Pure trap weakening corresponds to Vf < V,, while 
Lf = Li, and pure trap squeezing to Lf < Li, keeping 
Vf = Vi (Fig. [1]). We shall next describe the efficiency of 
atomic Fock state preparation by mixed trap reduction 
(Vf < Vi and Lf < Li) keeping constant the relative 
smoothness parameter at = at, and going from the Ui 
to the Uf families of traps. This procedure allows us to 
apply squeezing of the potential up to any desired value 
of Lf keeping its bathtub shape. 

Figure^iUustrates the full counting statistics of the re- 
sulting state in different limits of a trap reduction scheme. 
Both pure weakening and squeezing fail to produce an 
atom-number state since the condition in Eq. (O is not 
fulfilled. It was shown in [14[ that this limitation arises 
as the result of truncation of the final state both in coor- 
dinate (pure weakening) or momentum (pure squeezing) 
space, with respect to the desirable Fock state | JV = Cf). 
However, this state, whose full-distribution reduces to a 
Kronecker delta S n ,Cf (see Eq. (jTUJ) ) , can be obtained 
by combining both strategies, as shown in the middle 
panel. The final trap is then perfectly filled by the pure 
atom-number state. In what follows we shall characterize 
the efficiency of the method just by the mean and atom 
number variance of the prepared state, see Eqs. ([2]) and 
©. The normalized variance, crjj /{Nf) allows us to 
distinguish between the sub-Poissonian, Poissonian, and 
super-Poissonian statistics whenever it is lower, equal or 
greater than 1, respectively. 

Let us now consider different trap geometries. Gen- 
erally, the effect of the smoothness of the potential is 
to increase the density of states near the brims, where 
the spacing between adjacent energy states is reduced, 
see Fig. [3] As a consequence, a higher control of the 
depth of the potential would be required. Nonetheless 
Fig. H] shows that by increasing the smoothness, the 
Fock state creation condition is actually satisfied for 
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FIG. 3: Effect of the smoothness of the trapping potential 
(a) on the spectrum (b). As a increases the density of states 
concentrates near the brim of the trap and more bound states 
appear. The spectrum in the low panel is obtained for a 
potential V = (IOtt) 2 (in units of h 2 /2mL 2 ) . 
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FIG. 4: Asymptotic mean value (a) and variance (b) of the 
atom number distribution of a Tonks-Girardeau gas obtained 
by sudden weakening-squeezing as a function of the width 
ratio between the final and initial trap Lf /Li. The process 
is robust for different values of the smoothing parameter a. 
The initial state is assumed to be in the ground state. U a are 
chosen such that d = 100 and C/ = 10, this means that for 
the bathtub potential Ui = (IOOtt) 2 and Uf = (10tt) 2 , while 



for the Gaussian, Ui = (287r) 2 and Uf — 
we assume Ni — d = 100. 
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FIG. 5: Effect of temperature on the asymptotic mean value 
of the atom number distribution of a Tonks-Girardeau gas ob- 
tained by sudden weakening-squeezing for the case of a square 
potential as a function of the width ratio between the initial 
and final trap. In plot (a) Ui = (1007r) 2 and Uf = (107r) 2 so 
the capacities C a are the same as in Fig. [4] but the initial 
occupation is now chosen to be Ni = 0.8d. As we increase 
the temperature, fi/ksT < 5 the method starts to fail, but 
this can be improved increasing the capacity of the initial 
well, as shown in plot (b), where the weakening-squeezing 
process is applied at the same temperature but two differ- 
ent initial traps: Ui = (1007r) 2 used in the previous plot and 
U = (130tt) 2 , both with Ni = 0.8d. (Ei = h 2 n 2 /2mL; 2 .) 



We might conclude that an invariably efficient strategy 
for the sudden transition between a Ui and U / trap fami- 
lies, is achieved by reducing to half the width of the initial 
trap and reducing the depth to the trap accordingly, so 
as to achieve the desired Uf and capacity C/, 
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which warrants the preparation of the Fock state \Cf) 
corresponding to the [//-family. 



NON-ZERO TEMPERATURE 



a broader range of parameters which includes conditions 
nearer pure weakening and pure squeezing. This is be- 
cause the initial state is spread out along the same re- 
gion in configuration space as the final one; moreover the 
looser confinement reduces the momentum components 
of the final state, which can be resolved more easily by 
the initial state. 



The above formalism can be generalized in a straight- 
forward way to account for the atom number distribution 
resulting from an arbitrary initial state at non-zero tem- 
perature. It suffices to redefine 

A i = 5>„K)K| (16) 

n 
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which, in general, is not a projector now, where 7r„ is the 
occupation probability of the state \<p l n ). In the ground 
state of the TG gas n n = 1 Vn = 1, . . . , iV, and 7r„ = 
otherwise. For a thermal state, the Fermi-Dirac weights 
tt„ ={exp[^(£^-/i)] + l} _1 (with/3= l/fc B T where fc B 
is the Boltzmann constant and T the absolute tempera- 
ture) result due to the effective Pauli exclusion principle 
mimicked by bosons in the TG regime [23J. Notice the 
normalization ^2 n ir n = ^i- The preparation of a Fock 
state by a sudden change of the trap will still be feasi- 
ble as long as A/ C A t . The numerical results in Figure 
[5] (upper panel) illustrate the degradation of the quality 
of final state with increasing temperature. However, the 
lower panel shows that this negative effect of temperature 
can be compensated by starting from a "bigger" initial 
trap with larger capacity. 

VI. DISCUSSION AND CONCLUSION 

We conclude that the controlled preparation of atomic 
Fock states in the strongly interacting (Tonks-Girardeau) 
regime can be achieved by combining weakening and 
squeezing of the trapping potential. The process is robust 
with respect to the smoothness of the potential trap, and 
moreover the deteriorating effect of increasing tempera- 



ture can be compensated by enlarging the capacity of the 
initial trap. However, it is still an experimental challenge 
to get to the strong TG limit which must translate into 
a correction to the fidelity. By contrast, non- interacting 
polarized Fermions would be an ideal system for Fock 
state preparation. For ultracold fermions, due to the 
wavefunction antisymmetry, s-wave scattering is forbid- 
den and generally p-wave interactions can be neglected so 
that the gas is non-interacting to a good approximation. 
Such type of gases can be prepared in the laboratory 
with linear densities of the order 0.2 — 2 /im -1 for which 
the polarization remains constant in a given experiment 
[26| . For such gases the trap reduction technique can be 
directly extended to two and three dimensions. 



Acknowledgments 

We acknowledge discussions with Martin B. Plenio 
and Jiirgen Eschner. This work has been supported by 
Ministerio de Educacion y Ciencia (F1S2006-10268-C03- 
01, FIS2008-01236) and the Basque Country University 
(UPV-EHU, GIU07/40). The work of MGR was sup- 
ported by the R. A. Welch Foundation and the National 
Science Foundation. 



C. -S. Chuu, F. Schreck, T. P. Meyrath, J. L. Hanssen, G. 

N. Price, and M. G. Raizen, Phys. Rev. Lett. 95, 260403 [15 
(2005). [16 
I. Dotsenko et al, Phys.Rev. Lett. 95, 033002 (2005). 
N. Schlosser, G. Reymond, and P. Grangier, Phys. Rev. [17 
Lett. 89,023005 (2002). 

J. Sebby-Strabley, B. L. Brown, M. Anderlini, P. J. Lee, [18 
W. D. Phillips, and J. V. Porto, Phys. Rev. Lett. 98, 
200405 (2007). [19 
P. Bouyer and M. A. Kasevich, Phys. Rev. A 57, R1083 
(1997). [20 

D. Jaksch et al, Phys. Rev. Lett. 52, 1975 (1999) 

G. A. Prataviera, J. Zapata, and P. Meystre, Phys. Rev. 

A 62, 023605 (2000). [21 

J. Mompart et al, Phys. Rev. Lett. . 90, 147901 (2993) 

R. B. Diener, B. Wu, M. G. Raizen, and Q. Niu, Phys. [22 

Rev. Lett. 89, 070401 (2002). 

B. Mohring, M. Bienert, F. Haug, G. Morigi, W. P. Schle- 
ich, and M. G. Raizen, Phys. Rev. A 71, 053601 (2005). [23 
G. Nandi, A. Sizmann, J. Fortagh, C. Wei, and R. Walser, 
Phys. Rev. A 78, 013605 (2008). [24 
M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, I. [25 
Bloch, Nature 415, 39 (2002). 

A. M. Dudarev, M. G. Raizen, and Q. Niu, Phys. Rev. 

Lett. 98, 063001 (2007). [26] 

A. del Campo and J. G. Muga, Phys. Rev. A 78, 023412 



(2008). 

M. D. Girardeau, J. Math. Phys. 1, 516 (1960). 

D. S. Petrov, G. V. Shlyapnikov, and J. T. M. Walraven, 

Phys. Rev. Lett. 85, 3745 (2000). 

M. Olshanii, Phys. Rev. Lett. 81, 938, (1998); V. Dunjko, 

V. Lorent, and M. Olshanii, ibid 86, 5413 (2001). 

T. Cheon and T. Shigehara, Phys. Rev. Lett. 82, 2536 

(1999). 

M. D. Girardeau and E. M. Wright, Phys. Rev. Lett. 84, 
5691 (2000). 

L. S. Levitov and G. B. Lesovik, JETP Lett. 58, 230 
(1993); L. S. Levitov, H.-W. Lee, and G. B. Lesovik, J. 
Math. Phys. 37, 10 (1996). 

V. V. Kocharovsky, VI. V. Kocharovsky, and M. O. 

Scully, Phys. Rev. Lett. 84, 2306 (2000). 

I. Klich, in Quantum Noise in mesoscopic Physics, 

edited by Yu. v. Nazarov (Kluwer, Dordrecht, 2003); 

|arXiv:cond-"mat /0209642 

M. Budde and K. Molmer, Phys. Rev. A 70, 053618 
(2004). 

K. Schonhammer, Phys. Rev. B 75, 205329 (2007). 
See, for example, Numerical Recipes in FORTRAN: 
The Art of Scientific Computing (Cambridge University 
Press, 1992). 

M. Kohl et al, Phys. Rev. Lett. 95, 230401 (2005). 



